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ABSTRACT 



Context. M 1-92 can be considered an archetype of bipolar pre-planetary nebulae. It shows a clear axial symmetry, along with the 
kinematics and momentum excess characteristic of this class of envelopes around post-AGB stars. 

Aims. By taking advantage of the new extended configuration of the IRAM Plateau de Bure interferometer, we wanted to study the 
morphology and velocity field of the molecular gas better in this nebula, particularly in its central part. 

Methods. We performed sub-arcsecond resolution interferometric observations of the 7=2-1 rotational line of 13 CO in M 1-92. 
Results. We found that the equatorial component is a thin flat disk, which expands radially with a velocity proportional to the distance 
to the center. The kinetic age of this equatorial flow is very similar to that of the two lobes. The small widths and velocity dispersion 
in the gas forming the lobe walls confirm that the acceleration responsible for the nebular shape could not last more than 100-120 yr. 
Conclusions. The present kinematics of the molecular gas can be explained as the result of a single brief acceleration event, after 
which the nebula reached an expansion velocity field with axial symmetry. In view of the similarity to other objects, we speculate on 
the possibility that the whole nebula was formed as a result of a magneto-rotational explosion in a common-envelope system. 
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1. Introduction 

PN M 1-92, also known as Minkowski's Footprint, is an O-rich 
pre-pl anetary nebula (pPN) originally discovered by Minkowski 
( 1946), which has been deeply observed in the optical, NIR, and 
radio wavelengths. The central star has a spectrum compatible 
with T e ff = 6500 K and logg = 0.5, although at short wave- 
lengths it shows an excess that can be fitted with th e presence of 
a seco ndary star of T e ff = 18 000 K and log g — 5.0 ( Arriet aet al.l 
120051) . The distance to M 1-92 has been estimated as 2.5 kpc, 
adopting a nominal luminosity of 10 4 Lq for a post-AGB star 
(Cohen & Kuhill 19771) . The nebula has a size of 11" in the axial 
direction and 6" in the transversal one. It consists of a two-lobe 
reflection nebula, divided by a dense equatorial component. The 
two lobes define a clear axis of symmetry, oriented at a position 
angle (P. A.) of 311° (measured from north to east). This axis is 
inclined in respect to the plane of the sky by about 35°, with the 
northwest lobe pointing to us. Optical spectroscopy shows the 
presence of wide absorption components, revealing the existence 
of (bipolar) flows of ionized gas close to the star , with expansion 
velocities up to 750 kms -1 (lArrieta et al.l l2005). Line emission 
in Ha and in atomic forbidden lines ( O i, O in, N n, S ii) is de - 
tected from the middle of the two lobes (Bujarrab al et al.l!9 98b). 
probably tra cing the location of shocks p r opagat ing along these 
jets (see also lSolfll994T) . lBuiarrabal et al.l (fl998b) estimated that 
the mass of the ionized gas in the nebula is only about 10~ 3 Mq, 
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and lArrieta et al.l (|2P05) concluded that the present ionized mass 
loss is very low, 3 10 -13 Mq yr~' . 

Because of the early stage of the central star in its post-AGB 
evolution, almost all the nebular material is still in the form of 
molecular gas. This com ponent has been studied in detail by 
Bujarraba l et al.l (Qj?98a) by means of interferometric observa- 
tions of the 7=1-0 and 7=2-1 rotational lines of 13 CO, obtain- 
ing a resolution of ~1'.'0 at 220 GHz (in the 7=2-1 line). The 
molecular gas is located in the equatorial disk and in the walls 
of the two lobes, which are practically empty. From those maps, 
the total mass estimated for the molecular component is 0.9 Mq, 
0.2 Mq of which is located in the disk and 0.35 Mq in the walls 
of each lobe. The excitation was found to be very low across the 
nebula, from 10 to 15 K, and the densities ranging from 5 10 4 
to 3 10 5 cm 3 , the warmer and densest parts located close to the 
equator. The velocity field of the gas was also investigated using 
a spectral resolution of 3.25 km s~'. These authors found that 
the kinematics consists in a low radial expansion component 
of 8kms _1 , which dominates in the equatorial region, plus an 
axial component that increases linearly with the distance to the 
equatorial plane, reaching an expansion velocity of 70 km s~' at 
the two nebular tips. Such a Hubble-like velocity law was inter- 
preted as the result of a sudden interaction that lasted much less 
than the kinetic age of the axial outflow (~ 1000 yr). In addition, 
an axial momentum of 22 Mq km s was measured in the neb- 
ula. This value is much larger that what can be provided by radi- 
ation pressure on dust grains on the expected acceleration lime 
(about 100 yr or less). In fact these two properties, the presence 
of Hubble-like axial accelera tion, and the momen tum excess are 
very common among pPNe (Bujaii-ab ailtini200l . 
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Fig. 1. The 3.25 km s~' resolution channel maps of the intensity of 13 CO 7=2-1 in M 1-92 from the new+old Plateau de Bure data 
(see text). The LSR velocity in kms 1 for the center of each channel is indicated in the upper left corner. Drawn contours are -15 
(dashed), and 15 to 90 by 10 mJy beam 1 . The assumed coordinates for the reference position are indicated in the top left corner of 
the plot. 



In spite of all these observational efforts, very little is known 
about the mechanism that accelerated and stretched the nebula 
along its symmetry axis, as well as about what is powering the 
present bipolar outflow seen in ionized gas. (Note that the jet 
seen in the optical does not have enough energy or momentum 
to explain the present kinematics, so we may be in the presence 
of two different processes.) This is partially due to the limited 
spatial resolution of the molecular line observation, which can- 
not reveal the structure of the lobe walls or the shape and size 
of the equatorial component. Making use of the new, improved 
resolution of the Plateau de Bure interferometer, in this letter we 
report on new sub-arcsecond resolution observations of M 1-92 
in the 13 CO 7=2-1 line. 



2. New Plateau de Bure observations 

These new observations of M 1-92 were conducted with the 
new extended configurations 6Aq & 6Bq of the IRAM Plateau 
de Bure mm-wave interferometer. The observations were per- 
formed on January 24 and February 26, 2006 (configurations A 
and B, respectively), under excellent to good weather conditions. 
Using the 1.3 mm receiver, we observed the 7=2-1 rotational 
line of 13 CO, adopting a rest frequency of 230399 MHz. The 



7=1-0 lines of 13 CO and C ls O were simultaneously observed 
with the 3 mm receiver, but their results will be discussed else- 
where. The correlator was set up to cover the whole 7=2-1 line 
of 13 CO with an original spectral resolution of 156.25 kHz, i.e. 
0.20 km s _1 at 230.4 GHz; however, to ease the comparison with 
previous maps of the nebula, we re-sampled the channels to the 
same velocity resolution used before. To improve the reliability 
of the phase calibration, water vapor corrections based on the re- 
sults of the 22 GHz radiometers installed on the antennas were 
applied. Final phase calibration was performed from the obser- 
vations of the quasars J1923+210 and J2023+336 (inter-spaced 
every 20 min.), which were also used as secondary amplitude 
calibrators. Observations of the radio-sources 3C 273, 3C 345, 
and MWC 349 were performed to fix the absolute flux scale and 
to calibrate the gain across the band. 

After the data set was calibrated, maps were produced and 
CLEANed. The maps from the new data showed a significant 
fraction of lost flux, ~ 50% on average, when compared with 
previous observations. These losses are smaller, only ~ 30%, for 
velocities that are both extreme and close to the systemic one. 
In contrast, for intermediate expansion velocities, 5-20 kms~', 
the missing flux fraction is as high as ~ 2/3. In order to recover 
this lost emission and produce maps with both high-contrast 
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Fig. 2. Left: Position vs. velocity diagram for a cut along the symmetry axis of the nebula (P. A. 311°). Contours are as in FigQ] 
The straight line shows the velocity gradient found along this direction. Center: The same for the central region of the nebula (the 
rectangle shown on the left) from the maps obtained only with 6Aq data and for a spectral resolution of 0.5 km s -1 . Contours are 
as before. The white/black straight lines show the gradient of the equatorial/extended components of the nebula. Right: Selected 
channel velocity maps of the center of the nebula, obtained just from the 6Aq data with 0.5 km s -1 resolution. Contours are as before. 
The straight line delineate the axial direction. (See also Fig.Q]) 



and spatial resolution, we me rged our new data with those from 
our 1995-1997 observations (Buiarrab al et al.|[l998al) . using ro- 
bust weighting. The resulting clean-beam has an HPBW of 
0751x0735, with the major axis oriented at PA. = 36°. This 
clean beam was later used in the CLEANing procedure follow- 
ing the Clark method. Comparing the recovered flux in this data 
set with that from our previous observations both with Pico de 
Veleta and Plateau de Bure, we conclude that all spectra are 
compatible within the uncertainties and that no flux is lost in 
the maps. The final r.m.s. for the adopted spectral resolution of 
3.25 km s -1 is 4.5 mJy beam -1 , and the conversion from these 
flux units into Rayleigh-Jeans equivalent brightness temperature 
Tb is 142 K per Jy beam -1 . In addition we produced maps using 
just the visibilities from the most extended configuration (6Aq) 
and a spectral resolution of 0.5 km s -1 , to better investigate the 
structure and velocities in the equatorial component. The result- 
ing clean-beam has an HPBW of 0746x0726 at PA. of 25°, the 
final r.m.s. is 15 mJy beam -1 , and the conversion from these flux 
units into Tb is 208 K per Jy beam -1 . 

3. Observational results 

The maps are shown in Fig.Q] and the corresponding position 
vs. velocity diagram along the symmetry axis is displayed in 
Fig. |2] (left panel). The observations clearly show the structure 
of the nebula, with its two empty lobes divided by the equato- 
rial disk, and the axial Hubble-like velocity field. The lobe walls 
are very thin, with a de-convolved FWHP of 0725 (9 1 15 cm). 
The inner limits of these walls are much better defined than the 
outer boundary, where some weak diffuse emission is detected, 
especially at low velocities. The velocity dispersion in the walls 
must also be very small, since we can hardly detect it with a 
3.25 kins" 4 resolution. At the nebular tips, these widths and ve- 
locity dispersion are larger, about 079 and 5 km s -1 , respectively. 
All these figures are typically 10-15% of the distance to the cen- 



ter and of the expansion velocity at the corresponding points. 
This puts a strong upper limit on the duration of the forces re- 
sponsible for accelerating the nebula along the symmetry axis; 
we confirm that the process could not last more than about one 
tenth of the kinetic age of the flow: i.e. £ 100-120 yr (see later). 

At the center of the nebula we find that the equatorial struc- 
ture is not compatible with our previous model. We do not see 
the signa tures of a disk e xpanding at a constant radial velocity 
(see e.g. IChiu et alj|2006h . This is much better displayed with 
the results from the 6Aq configuration (Fig.|2j center and right 
panels). In the axial velocity cut, we find that this structure also 
shows a constant velocity gradient, but in the opposite sense. 
In the maps, we see that the emission at each velocity traces a 
strip parallel to the symmetry axis. All these features can only 
be fitted with a disk in which the radial expansion is not constant 
but proportional to the distance to the center: i.e. we also find a 
Hubble-like velocity field in the equatorial plane, a totally unex- 
pected result! Note that the maximum in the maps arises from 
the center of this disk and at the systemic velocity, in total dis- 
agreement with the expectations for an equatorial undisturbed 
remnant of an AGB envelope. From the extent of the emission at 
the systemic velocity, we derive that the diameter of the disk is 
174 (5 10 16 cm) and the height at half power is < 072 (7 10 15 cm). 
At the resolution achieved in the maps, we do not find a central 
cavity or signatures of gas rotation. 

For the axial and equatorial outflows, we find a projected 
velocity gradient of 7.6 and 12kms -1 per arc-second, respec- 
tively. Assuming the distance and inclination with respect to the 
plane of the sky given in Sect. 1, these figures translate into 
kinetic ages, t^ m , of 1060 and 1370yr, respectively. However, 
since the two flows are perpendicular to each other, these kinetic 
ages have an opposite dependence with the adopted inclination 
angle. In particular, changing this angle from 35° to 38.5° (i.e. 
only by 10% and within the uncertainties in our previous mod- 
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eling), both ages become the same, 120 yr. Th is inclination is 
still compatible the one derived bv lSolfl d 19941) of 33° ±5°, and 
the resulting coeval scenario is much more plausible than hav- 
ing one ejection along the equatorial plane 300 yr before the the 
axial flow is launched. If this hypothesis is correct, the whole 
CO nebula would share the same Hubble-like velocity field: the 
velocity would just be proportional to the distance, V = R/t^, 
and we can simply use the velocity along the line of sight to 
derive the gas distribution in this direction. In this case the po- 
sition vs. velocity digram shown in Fig. [2] would represent the 
spatial structure of M 1-92, if we forget our finite resolution in 
velocity and space, by simply translating the angular distance 
and velocities into length units using the distance to the target 
and the derived velocity gradient. Note that even the location of 
the weak diffuse emission seen along the equatorial plane at ve- 
locities around 10-15 km s _1 (Fig.|2]left) seems consistent with 
the halo detected in the maps in Fig.Q] 

4. The origin of the pPN M 1-92 

The origin of bipolar outflows in pPNe is not yet understood. 
Current scenarios involve the presence of accretion on the cen- 
tral star or on the companion in the case of binary systems 
through a rotating-disk, mag neto-centrifugal launc hing, and a 
variety of collimating agents (Balick & Frank 2002). The ques- 
tion of whether close binary systems are mandatory for pro- 
ducing axis-symmetric nebulae is also a long-standing issue 
(BujarrabaTetaDHoOO). Of course explaining the formation of 
a nebula involving coeval axial and equatorial Hubble-like flows 
becomes even more difficult. 

However, M 1-92 is not the only object in which complex 
structures resulting from a single sudden event have been found. 
The brightest luminous blue variable, r\ Car, is surrounded by a 
cloud of gas and dust, whose central and densest part is known as 
the "Homunculus" . This nebula, with an estimated mass ~ IMq, 
consis ts of a two- l obe s tructure divided by a thin equatorial 
torus. iMorse et"aH (1200 ll) performed a proper motion study in 
this nebula by comparing HST images taken at three different 
epochs and concluded that the "Homunculus " also shows a sin- 
gle Hubble law and was formed in less that 10 yr during the 
great eruption in the 1840's. These measurements are only done 
in the plane of the sky (neither the position nor the velocities 
are derived along the line of sight, as in our case), so they do 
not depend on the orientation of the nebula or on the distance. 
Although T] Car and the central system in M 1-92 are very dif- 
ferent, the striking similarity found in the kinematics and shape 
of these two ne bulae strongly su ggests that they may share a 
common origin. iMatt et all ((2006) have developed a model for 
j] Car in which the differential rotation between the stellar nu- 
cleus and its envelope can drive a magneto-rotational explosion. 
In this model, as a result of the sudden release of the magnetic 
energy when it overcomes the weight of the upper layers, a bipo- 
lar jet is launched but an equatorial flow is driven too. These au- 
thors suggest that this mechanism would also explain the shape s 
of bipolar pPNe like IRAS 1710 6 -3046 (IMatt & Bahddl2004h . 
Recently iNordhaus & Blackmanl (120061) conclude that for this 
mechanism to operate, how the differential rotation is sustained 
against the increasing magnetic force can be explained by the 
angular momentum transferred to the stellar envelope by a low- 
mass companion during a common envelope phase. If this is 
also the case for M 1-92, we should conclude that this pPNe 
is largely formed by gas directly expelled from the stellar enve- 
lope and not a former AGB circumstellar envelope accelerated 
by the interaction with some post- AGB flow. Otherwise it would 



be very difficult to explain how the present gas could show such 
an accurate ballistic movement. 

Although M 1-92 is so far the only case in which we can 
argue that the central disk was formed at the same time as the 
lobes, there are other pPNe in which radial acceleration has been 
found along th e equator i n ad dition to the typical bipolar flows. 
In CRL 2688 ICox et al.1 (l2000h found in the CO maps up to 7 
different bipolar outflows, three along the axis of symmetry (D 
to F in their naming) and four more (A to C and G) that could be 
running in different directions in the equatorial plane. In all these 
jets the velocity increases with the distance to the center. The au- 
thors do not comment on whether a Hubble law also holds here, 
but from Fig. 3 in that paper we conclude that a constant veloc- 
ity gradient is consistent with the data. However, it is very dif- 
ficult to see whether these equatorial flows have similar kinetic 
ages (so as to compare them with those for the axial jets). For 
the axial flows , this Hubb l e-law velocity field has been recently 
confirmed by lUeta et all d2006l) . by means of proper motion 
measurements using H ST NICMOS images. Another i nteresting 
case is CRL 618. Here Sanc hez Contreras et al.l ([2004 ) observed 
I2 CO, and when fitting their results with a numerical model, they 
found that an increasing velocity was also necessary for the cen- 
tral dense core. This model neither assumed an exact Hubble law 
field nor constrained this kinematics to a thin equatorial disk, but 
note that the fitting was not very precise in this region because of 
the limited spatial resolution of the observations. (In these two 
cases, the interaction of the post-AGB flows with the old AGB 
shell is clearly detected.) Finally, Frosty Leo presents a dense 
equatorial disk/torus that expands at higher velocities than those 
typical in AGB envelopes. Castro-Carri zo et al.l d2005i) conclude 
that such a high velocity was probably due to interactions be- 
tween the AGB and post-AGB flows in the a plane perpendicular 
to the axial direction. However, in view of our current finding, 
we may also speculate with the hypothesis of an equatorial ac- 
celeration in this source, too. Note, however, that M 1-92 is one 
the very few pPN in which the molecular gas has been studied 
at such high spatial and spectral resolution. Therefore it would 
very premature to conclude whether the mechanism proposed 
here for explaining the kinematics of M 1-92 should also apply 
to other pPNe. We must wait until high-resolution images of the 
molecular gas in more sources are available, but we think that 
this possibility should be explored in the context of understand- 
ing the issues of pPNe and PNe formation and shaping. 
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